arXiv: 1501.0459Ivl [gr-qc] 19 Jan 2015 


Rapidly rotating neutron stars in R-squared gravity 
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f{R) theories of gravity are one of the most popular alternative explanations for dark energy and 
therefore studying the possible astrophysical implications of these theories is an important task. In 
the present paper we make a substantial advance in this direction by considering rapidly rotating 
neutron stars in R^ gravity. The results are obtained numerically and the method we use is non- 
perturbative and self-consistent. The neutron star properties, such as mass, radius and moment of 
inertia, are studied in detail and the results show that rotation magnifies the deviations from general 
relativity and the maximum mass and moment of inertia can reach very high values. This observation 
is similar to previous studies of rapidly rotating neutron stars in other alternative theories of gravity, 
such as the scalar-tensor theories, and it can potentially lead to strong astrophysical manifestations. 

PACS numbers: 


I. INTRODUCTION 

The interest in studying the modifications of general relativity is growing significantly in the past few years. From 
one side this is connected with our aim to better understand the theory of gravity and to explore what are the possible 
deviations. But on the other hand there are observations which do not completely fit into the standard framework 
of Einstein's theory of gravity, such as the dark energy phenomena. Exactly this phenomena is the main motivation 
for the recent interest and advance of /(R) theories of gravity which can offer an alternative and viable explanation 
for the dark energy. 

Naturally/(R) theories of gravity were mainly investigated in a cosmological context. But every theory of gravity 
has to be able to also pass various astrophysical tests. As far as black holes are concerned, no-hair theorems exist 
M, according to which the solutions in /(R) theories are the same as the ones in general relativity (in the case of 
linear matter sources). But this is not the case with neutron stars, where the presence of matter can lead to significant 
deviations. On the other hand there is a variety of observations of compact stars that can be used to test both the 
weak field and the strong field regime of different theories of gravity. That is why studying neutron stars in /(R) 
theories and confronting them against the observations, is an inseparable part of the attempts to find a generalized 
theory of gravity that can naturally incorporate the dark energy phenomena. 

Compact stars in /(R) theories were intensively studied in the past few years HQ. Due to the complexity of the 
field equations, the more detailed investigation of theproblem and the confrontation with observations were done 
primarily using a perturbative approach (see e.g. |3-0)- But it was later shown that studying the neutron stars in 
f(R) th eories of gravity perturbatively is misleading and one has to solve the full field equations self-consistently 
Ill4l.llal . As a matter of fact the calculations can be simplified considerably if one uses the mathematical equivalence 
between the /(R) gravity and a particular class of scalar-tensor theories with a nonzero potential for the scalar field. 

The examination of nonrotating compact stars in /(R) theories showed that constraints on the parameters can not 
be set via the present observations of the neutron star masses and radii alone, since the uncertainties in the nuclear 
matter equation of state can be bigger than the deviations from the pure general relativity llT3l . But the results in 
the case of slow rotation suggest a possible test of /(R) gravity via the expected future observations of the neutron 
star moment of inertia that can impose tight constrains on the free parameters in the theory fl^ . A natural 

extension of these results is to consider rapid rotations. An astrophysical motivation comes from the observations of 
millisecond pulsars and the models of supramassive neutron stars formed after a merger, where the rotational effects 
will be non-negligible. On the other hand it was shown in lIT^ that rotation can significantly increase the deviations 
from general relativity in the case of scalar-tensor theory. Taking into account that /(R) theories are mathematically 
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equivalent to them, we can expect a similar magnification of the differences, thaf can potentially lead to additional 
observations constraints on f{R) theories at astrophysical scales. 


II. BASIC EQUATIONS 

We will describe the analytical setup of the problem very briefly. The reader is referred to for a detailed 

discussion of the problem. 

The action in f {R) theories of gravity is given by 

S = j d‘^Xy/^f{R) + Smatter{giw,x), (1) 

where R is the Ricci scalar curvature with respect to the space-time metric gf,v, Smatter is the action of the matter, and 
the matter fields are collectively denoted by x- The inequalities <f-f/dR^ > 0 and df / dR > 0 have to be imposed in 
order for the theory to be free of tachyonic instabilities and the appearance of ghosts. In the present paper we will 
be concentrated on the so-called R?- theory of gravity for which 

f{R)=R + aR^, (2) 

where a is a paramefer satisfying a > 0 in agreement with the above mentioned inequalities. 

In our calculations we will use the fact that the f{R) theories are mathematically equivalent to the Brans-Dicke 
theory with a parameter cobd = 0 and nonzero potential of the scalar field. The Jordan frame action is given by 

S = J ~g ~ bl(<l>)] -|- Smatter(^^i//i^)/ (3) 

where the gravitational scalar <3> and the potential 11(0) are defined by O = and 11(0) = — f{R), respec¬ 

tively. In the particular case of R^ gravity, that we are considering, one can easily show that <3> = 1 -|- 2aR and the 
scalar-field potential is Lf(0) = ~ 

Even though the Jordan frame is the physical frame where the physical quantities are measured, in practice it is 
often more convenient to study the scalar-tensor theories in the so-called Einstein frame, where the metric is 
defined by fhe following conformal transformation = Og^v The action in Einstein frame can be written in the 
form 


s = J d^x^* [R* - 2g*’^%cpd,cp - V{cp)] + S^atter(e (4) 

where R* is the Ricci scalar curvature with respect to the Einstein frame metric g))^ and the new scalar field cp is 

defined by ^ InO. The Einstein frame potential V{cp) is correspondingly V[(p) = A‘^{cp)U{<i>{(p)) where the 

Einstein frame coupling function A{(p) is given by 

a2(^) = (5) 


Eor the particular case of R^ gravify the scalar-field potential takes the form 


V(cp) 


1 

4fl 



-Al 
e V3 


2 


( 6 ) 


Due to the significant simplification of the field equations in the Einstein frame we will employee it^. Since this is 
not the physical frame, one has to transform the quantities, such as mass, radius and moment of inertia, back to the 
physical Jordan frame. A more detailed discussion of the problem can be found in fT3.IT^[T§l . 


^ The Einstein frame quantities will be marked with an asterisk. 
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FIG. 1: The mass of radius relation for EOS APR4 (left panel) and SLy4 (right panel) in the case of static neutron stars and stars 
rotating at the Kepler limit. Different styles and colors of fhe curves correspond to different values of fhe parameter a. 


Since we will consider rotating neutron stars, the following general ansatz for the stationary and axisymmetric 
Einstein frame metric can be used 

dsl = sin^ 6{d^ — codt)^ + e^‘^{dr^ + r^dO^), (7) 

where all the metric functions depend on r and 6 only. The dimensionally reduced field equations for the metric 
functions, the scalar field cp and the equation for hydrostationary equilibrium can be found in the Apendix (see 
also 13). The quantities, such as mass M and angular momentum /, are calculated using integrals throughout the 
volume of the star and the details are given in M- The moment of inertia of the star is given by the standard formula 
I = //n, where Q is the angular frequency. 

The field equations are solved using an extended version of the RNS code ll^ where the modifications coming 
from the scalar-tensor theory are implemented 13 . But the code developed in Il3 is only for the case of zero scalar 
field potential. Here we extend it by taking into account the nontrivial scalar field potential V{(p) which is directly 
connected to the choice of the function f{R). Such an extension is nontrivial, because the introduction of a potential 
makes the scalar field equation very stiff. That is why one has to be very careful with the choice of the auxiliary 
parameters, such as the numerical accuracy and the relaxation factor, in order to be able to converge to the desired 
solutions. In addition to the tests performed in 13 . we have verified our code against the results obtained for static 
and slowly rotating neutron stars in gravity llT3.f3 . 

In the next section, where the numerical results are presented, we shall use the dimensionless parameter a —> a / Rg 
and the dimensionless moment of inertia I —)• I/MqRq, where M© is the solar mass and Rg is one half of the solar 
gravitational radius Rg = 1.47664 km (i.e. the solar mass in geometrical units). 


III. NUMERICAL RESULTS 

In our studies we have chosen two realistic equations of state (EOS) that fulfill the current observational constraints 
on the neutron star masses and radii i20| - l3 . These are the APR4 ll2^ and SLy4 EOS ll3l . This choice is motivated 
also by the fact that both of them were used in our previous studies of static and slowly rotating neutron stars in 
/(R) gravity llT3.[3 . that eases the comparison of the results. 

In Eig. [T]the mass of radius relation is plotted for several values of the parameter a in the R^ theory of gravity. The 
limit fl —>• 0 corresponds to pure GR and the case with a = 10^ represents nearly the maximum possible deviation 
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FIG. 2: The moment of inertia as a function of mass for EOS APR4 (left panel) and SLy4 (right panel) in the case of static neutron 
stars and stars rotating at the Kepler limit. Different styles and colors of the curves correspond to different values of the parameter 
a. 


from GR (values of a > 10^ give resulfs very close to the a = 10^ case). Two limiting sequences in terms of rotation are 
shown - the nonrotating neutron star sequences and the case of stars rotating at their Kepler limit, which represents 
the maximum possible deviation from the static solutions. As can be seen on the figure, the rapid rotation magnifies 
the differences between f{R) theories of gravity and GR - the maximum masses increase considerably for large a at 
the Kepler limit. The situation is similar to the case of scalar-tensor theories of gravity with zero potential for the 
scalar field which admit scalarization of the solutions . In these theories considerable differences with GR are 
observed for rapid rotation while in the static case the scalarized solutions are very close to GR. In addition, the range 
of parameters where such scalarization exists is significantly broadened for fast rotation. Since the /(R) theories of 
gravity are mathematically equivalent to a particular class of scalar-tensor theories, similar increase of the deviations 
compared to the static limit is expected and justified. 

As the studies in the slow rotation regime showed, the gravity has more pronounced effect on the rotational 
properties of the star such as the neutron star moment of inertia jl^ . The differences with GR can even exceed the 
EOS uncertainty for large values of the parameter a that can be used in order to impose observational constraints on 
the theory. As expected, the rapid rotation magnifies even further these deviations, similar to the case of the M(R} 
dependence. Indeed this can be seen on Fig. |2]where the moment of inertia is plotted as a function of mass for static 
model and for sequences of models rotating at their Kepler limit. 

A quantitative estimation of the relative deviations of the maximum mass and the maximum moment of inertia 
are given in Table I for the two EOS we are considering. The differences are given in percentages and for a fixed 
value of a the quantities AMmax and Atmax are defined in the standard way 


AAImax 


^«max -'"max 


AIrnax — 


t/(R) _ tGR 
imax ^max 

jGR 
^max 


( 8 ) 


Our results show that AMmax and Atmax are very small for values of a < 1. As a matter of fact if one continues 
decreasing the parameter a, these deviations can become even negative as shown in (13.IT^ . Therefore the maximum 
mass in R^ gravity can drop below the GR one for a < 1 (for the EOS we are considering), but the deviations 
AMmax and Atmax will have very small absolute values below 1%. Such differences can be hardly used to set any 
observational constrains on the f{R) theories and that is why we would not pay special attention to the case when 
a < 1. 

As the results show, the maximum mass increases up to 10% in the static limit and above 16% in the Kepler limit. 
As already commented in Q such deviations are non-negligible but they are comparable with the uncertainties in 
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the nuclear matter equation of state. That is why it would be difficult to set constraints on the f{R) theories using 
measurement of the neutron star mass and radius alone, until the EOS can be determined with smaller uncertainty. 
But one should keep in mind the following. Some of the proposed methods for consfraining the nuclear matter 
EOS are related to indirect measurement of the maximum mass of rapidly rotating neutron stars. More precisely, 
as the results in |2^ show, the observation or the lack of observation of gravitational waves signal emitted by 
oscillating post-merger neutron stars (a few milliseconds after the merger) can be used to determine with a good 
accuracy the upper limit of the mass of stable rapidly rotating neutron stars. But exactly this is the regime where 
the f{R) gravity will give the most significant deviations. That is why all strategies for determining the nuclear 
matter EOS should be used with caution since modifications of fhe gravitational theory can produces similar (both 
qualitatively and quantitatively) effects, especially in the rapidly rotating regime where the deviations from GR are 
substantially magnified. 

The differences in the neutron star moment of inertia on the other hand can be much more dramatic - from 
40% increase of the maximum angular momentum in the static limit, to over 65% in the case of rapid rotation. As 
commented in | [l^ such lar ge d eviations can be potentially measured by the forthcoming observations of the neutron 
star moment of inertia ll^ that can lead to a direct test of the R^ gravity. The gravitational wave spectrum of 
rotating neutron stars would be also altered considerably and such a study is underway. In particular the criteria 
for the development of the secular Chandrasekhar-Friedman-Schutz instability, what sets in due to the emission 
of gravitafional radiation, can be altered significantly that would have a direct effect on the predictions for the 
detectability of the signal. The validity of various phenomenological EOS independenf relations, connecting the 
neutron star mass, radius, moment of inertia, etc. 1^1^ , including the I-Love-Q relations IMnla]/ has to be also 
checked and such a study in underway. 


IV. CONCLUSION 

In the present paper we have studied rapidly rotating neutron stars in R^ theories of gravity by throating the 
problem non-perturbatively and self-consistently. For this purpose we have used the well known fact that f{R) 
theories are mathematically equivalent to a particular class of scalar-tensor theories with nonzero potential for the 
scalar field. The numerical solutions are obtained using an extended version of the RNS code that can deal with such 
scalar-tensor theories. The equilibrium properties of the solutions are studied and a detailed comparison with the 
static case is performed. 

The results show that similar to the neutron stars in scalar-tensor theories, the rotation can magnify significantly 
the deviations between the f{R) neutron stars and the corresponding solutions in pure GR. This effect is better pro¬ 
nounced for fhe rotational properties of the star, such as the moment of inertia and the angular momentum. The 
differences can be comparable or larger fo the rmcertainties in the nuclear matter EOS. The maximum mass of neu- 
fron sfars in R^ gravity can increase more than 16% in the case of rapid rotation, while the increase of fhe maximum 
moment of inertia reaches 65%. Clearly such large differences can have a non-negligible influence on various astro- 
physical phenomena, such as the observed quasiperiodic oscillations, the neutron star oscillation spectrum and the 
associated gravitational wave emission, etc. The deviations are in general larger for larger values of the parameter 


TABLE I: The relative deviations of the maximum mass and the maximum moment of inertia. 



APR4 Static 


APR4 Kepler 


a 

AAlnTax[%] 

AImax[%] 

AMmax[%] AJnTax[%] 

GR 

0.0 

0.0 

0.0 

0.0 

1 

0.1 

2.9 

1.1 

7.1 

10 

2.6 

15.6 

6.3 

28.1 

10^ 

5.7 

28.5 

11.5 

50.2 

10^ 

7.6 

36.9 

14.4 

63.8 


SLy4 Static 


SLy4 Kepler 


a 

AMn.ax[%] 

AImax[%] 

AMmax[%] AW[%] 

GR 

0.0 

0.0 

0.0 

0.0 

1 

0.9 

4.6 

2.1 

7.8 

10 

4.5 

18.4 

8.0 

29.6 

10^ 

7.8 

32.1 

13.5 

52.4 

10^ 

10.0 

41.2 

16.6 

65.6 
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a in the gravity, that can be used to impose an upper limit on this parameter. The various EOS independent 
relations between the neutron star properties might be also modified, carrying a signature from the modified theory 
of gravity. Such a study is underway. 

Once the high density equation of state is better know, the results in the present paper can se used to set tight 
constraints on the f{R) theories. But since rapid rotation can be observed for various astrophysical objects, such as 
the millisecond pulsars and the supramassive neutron stars formed after a merger, various strategies for extracting 
the neutron star EOS should be applied with caution. The reason is that modified theories of gravity can lead to 
competitive effects compared to a change in the EOS, especially in the rapidly rotating case. 


Acknowledgments 

D. D. would like to thank the Alexander von Humboldt Eoundation for support. K.K. and S.Y. would like to thank 
the Research Group Linkage Programme of the Alexander von Humboldt Eoundation for the support and S.Y. would 
like to thank the Institute for Theoretical Astrophysics Tuebingen for its kind hospitality. The networking support 
by the COST Action MP1304 is gratefully acknowledged. 


Appendix A: Dimensionally reduced field equations 

In this appendix we give the Einstein frame dimensionally reduced stationary and axisymmetric field equations 
describing the structure of the rapidly rotating neutron stars in f{R) theories. The Einstein frame fluid energy 
density, pressure and 4-velocity will be denoted by e*, p* and m* . The same quantities in the Jordan frame will be 
denoted by e, p and uf‘. The relations between the quantities in both frames are explicitly given by e* = A^(^)e, 
p* = A^{(p)p and = A{cp)ut‘. 

In writing the dimensionally reduced equations it is convenient to use the proper velocity v of the fluid given by 


V = {Ci —cv)r sin 6e (Al) 

where O is the fluid angular velocity defined by O = -^. In the present paper we consider only uniformly rotating 
neutron stars with constant angular velocity. O and v are the same in both Einstein and Jordan frame lIT^ . The fluid 
four velocity in the Einstein frame then is 

,, e-(v+T)/2 

ul= , -- [1,0,0,0]. (A2) 

V 1 — 

Erom a numerical point of view it is more convenient to use the angular coordinate }i = cos 6 instead of 6. In terms 
of this coordinate, the dimensionally reduced field equations for the metric functions 7 , cr and co are the following 


(a + b, - P 


£7/2 |(1577-p^ _ yj(pjje2a_|_ 



(Ib/rp* 




1 

2 




1 1 - p 2 

2 r 2 



(A3) 


A{cre'^^^) = js/rje*-h -h 7^(1 — p^)e 


-2cr 


{drU))^ + ^ 


I. P 3 


a 



(Ib/rp* 


V{cp))e^^ 


-drj- 


F . ^2 ^2 


(A4) 
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1 . A 


(e* + p*)(n-a;)^ 2 a , 

1 - i ;2 + 

F. .1 


1 l-p 2 


- -a. (-7 + 2^) + (-7 + 2c.) - - (a.7)^ - 4 ^ (3,7)^+ 

+ (^rcrf + - ''^(1 - (^( aro ;)^ + ( 3 ^ 0 ;)^^ 


_ 8 ^ £.(l + ^)+ 2 p,t;^ ^ 2 . _ 1 y(<p), 2 « 

1 - u2 2 


Here, the (differential operator A is defined by 


, ^2 1 1 “ -^2 2 m 

^ ~ + 7 ^'' ): 2 — 


(A5) 


(A6) 


For the metric function a we have two first order partial differential equations. In the numerical method in order 
to determine a we need only one of them, namely the following one 


a^a = - 


3^7 + ^//^ - {(1 - p 2 )(l + J ^ 


(A7) 


2 1 


rdrirdy'Y) + r^{Mf - (1 - F^)(3f/7)^ - - F^)a^7] + }id^j x [-p + (1 - p^)a^ 7 ] + 


+ 7[“F + (1 - X l^idrj+ dr(T)^ - (1 - + d^yCrf + AA{d(p)^ - 4(1 - jA){djy(pf 


+pra,- 7 [l + rdyj] — (1 — + ^dyj) 


diidyj + 3 ^ 73^7 + + djia){dy^ + d^a) + IdjyCpdyCp 


2 

+ ^(1 - [-[-}i + (1 - }i^)d^j][A{drCv)^ - . 2(1 - ii^){d^cv)^] + 

+2{1 - }A)AdnLodrCo{l + rdrj) 

The dimensionally reduced equation for the scalar field cp is 


+ 


+ 


1 — m2 

Acp = -drjdyCp - —d^jd^(p ■ 


V3^* dcp J 


.,2a 


The above equations have to be supplemented with the equation for fluid equilibrium 


(A8) 


^iP 
e + p 


a,(In ui) + 


= 0 . 


(A9) 
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